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Abstract 

Analogies between Coulomb excitations of nuclei and ionic molecules 
by charged projectiles is utilized to calculate vibro- rotational 
excitations of E 2 h molecular ions by e + impact by a semi- classical method 


developed in the nuclear case 
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1 . Introduction 

Recently there has been a considerable upsurge of interest in the study 
of rotational and vibrational energy loss processes in diatomic molecules 
by electron impact. They are of much importance not only for understanding 
the fundamental energy exchange processes involved, but also for their 
applications in such allied fields as Astrophysics and Atmospheric Physics. 

In the present work we shall investigate the coupled excitation of 
vibro-rotational states of hydrogen molecular ions, Hg + , by collision with 
positrons, e + . The study of such excitations with e + is not only important 
for its intrinsic significance but also for the mathematical simplicity 
it introduces in the formulation of the complex excitation process itself. 

This is due to the fact that the Pauli exchange does not enter directly 
into the problem. 

In this work we shall adopt a semi-classical view and make use of the 
analogy of Coulomb excitations of nuclei, which has been studied extensively 
in the past (Alder et al., 1965)* The present method is semi- classical 
in that we shall treat the target system quantum mechanically while the motion 
of the projectile would be assumed to be along a classical Coulomb trajectory. 
In the end, however, we shall attempt to modify the classical nature of the 
projectile motion, which does not distinguish between the initial and final 
states, by demanding that the principle of lurdtyoocbiybe satisfied by the 
cross-sectional, expressions and invoking the Correspondence Principle to be 
applied to such quantities as the classical velocity of the projectile. 
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2. Theory 

Let N, n and J denote, respectively, the electric, the vibrational and 
the rotational quantum numbers of the target molecule. In this paper we 
shall confine ourselves to the ground electronic state of the target and 
mostly disregard the quantum number N. One of the basic assumptions of 
our semiclassical approximation is that in the first approximation the 
incident positron moves essentially along a Rutherford trajectory around 
the effective positive charge of H 2 + , situated at the qp. of the molecule. 
Thus, while the positron moves along its trajectory it is allowed to induce 
transitions in the molecular motion through the electric coupling with the 
various molecular charges. The transition probability, b^ f , for the initial 
state i going to the final state f may then be given by the first order 
time dependent scattering theory of Dirac (I 9 B 6 ). Thus we write 




( 1 ) 


where iw = Ei - is the energy difference between the states i and f and 
Vint (t) is the interaction potential between the incident. positron and the 
H 2 + target. It can easily be seen (Fig. 1) that 
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where r = r(t) is the trajectory of the positron and rx is the position 
vector of the target electron, measured from the c.m. of the molecule. The 
vector R stands for the separation between the two nuclei of the target. 

We shall describe the target molecule in terms of simple product of 
normalized wave-functions (r x ), X^(R) and Yj m J (r) corresponding to the 
electronic, the vibrational and the rotational motions. Thus the total 
target wave function becomes 


|N„j>r $ («;>*„(*) Yj"(S) 


(3) 


In view of the fact that the incident positron repels itself from the 
positvely charged target, we may for sufficiently low energies, simplify 
the potential interaction (2) by expanding it for essentially non- penetrating 


projectile orbits. Thus we find 

A = *A 

Vt'n-t (*) - 


Y*x « ) Y: (i) 

l .,1, ^ 

! Y* (8) W ] 


A = - * 
As l,4V 

A =" X 


(J £A) 

h*" 


K >K. 


Taking the matrix elements between the initial state |N^ > and 

the final state |N f n f j f > we find from (3) and (4), 


A 1 *a 


(4) 


< 4 i vw(*>n> = { - ^ \ 


w 
A =-A 


( 5 ) 
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-* o£c* > h£(fti'C(^ x ] 

We find that in the present approximation, the above equation(5) controlls • 
the entire process of transitions among the various states and the various 
modes of motion of the target molecule. It can easily be seen that the 
perturbation of the electronic motion, represented by the first sum in (5) 
is independent of the nuclear perturbation, represented by the second sum 
in (5). This is, of course, a consequence of the BornjOppenhelmer separation 
already assumed in the product wave function . (5). We note however that the 
transitions among any two vibrational states ln^ > and n^> (n^ / n^) or 
any two rotational states |j^> and |j f > (J^ ^ J f ) or both lead to a vanishing 
of the first sum in ( 5 ). Thus the vibrational- rotational transition probab- 
ilities in the ground electronic state is found to be given by 

/in X 

\>u - - 4ir«; 2. 

/L* - * 

= \ e l lit \K y* ( *<%>; <k > 


nj> = 


5 e‘ [k «)}-*■' K ( <(«; ck , 
I X„ («) Y, 0 ’ ( &> 


We note that all informations regarding the trajectory of the incident 
position is contained in the orbit integrals J\ IL • Evaluation of these 
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intsgrals is extensively discussed in the literature (e.g. Alder et al., 1965) 
and we shall note here that the Coulomb trajectory of e"*" can most conveniently 
be given by the parametric equations 


r(t) = r Q (I cosh t + 5.) 

x(t) = r (cosh t + E) 
o 

y(t) = r (E a - l)^ einh t 
o 

z(t) * 0 

t * — 0 (E sinh t + t ) 
u o 


where E is the eccentricity, t is the eccentric anomaly, r Q is half the 


distance of closest approach, u Q is the projectile velocity and t is the 
time parameter. Substituting (7) In ( 6 ) and choosing the quantization 
axis along the angular momentum of the molecular rotation, it can be 
shown (Alder et al., 1965) 

X Y * ^ T/ (T,o) 


Ix* * 5 e 


i * 1 (t *t Q+ r v )^ siwU x 




( f + CbahT + <: (rf-O 1 x ) 

( I 4- XCbshT ) ^ 


dx 


( 8 ) 
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co c = 


!=l - f * 


60 , 


Po 

V-O ) 
I 

2 -f= 


( 9 ) 
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We have defined E Q to be the incident energy in a.u. and 6 as the 
scattering angle. To calculate the cross-sections we need the number of incident 
particles in a plane wave of unit flux with impact parameters between b 
and b + db. This is given by 

jutolL - — — x d-H. where df] is an elementary 

■«(*■’*' solid angle ; 

k zs K* Gsfc 


The differential cross-section for the \ multipole transition averaged 
over the initial substates m^ and summed over the final substates m^, may 
therefore be written as 


©»• n L -.,5UWl‘ 


■J-kM . 4rl ^ [ r>* (;,■)!»/. 

d_a ^2> + )) 4 ' Ar ‘ x ^4)4 


The total cross-section is obtained by integrating (11) over all scattering 
angles. Thus 


- i,'.U 
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In (11) and (13) the quantity is the reduced transition probability which 


we have defined to be 


BxCnJf-fcM* )= 21 |<^+l W A Y*(f)| 

A ** 


I Yx lh^> I 


where 


<**h ii<?/otfku> 

= a; x A x (n,^n 4 )<iJ y^||S c > 


with the vibrational matrix elements 


ft nine •*",)« ^ \(^ ** VM* 


and the reduced matrix element 


<s<41 Y*ll 1> = (- 1 ) 1 


(de-Shalit and Talmi, 1963) . 
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3. The Symraetrization of Classical Cross-Sections 

In the present approach, the projectile orbit has been described 
classically and consequently the cross-sectional expressions ( 11 ) and ( 13 ) 
do not distinguish between the initial and final states of the trajectories. 

One of the most successful ways of introducing the initial and final parameters 
for the orbits is to impose the principle of reciprocity on the crose- 
sectional expressions and replace various classical parameter? by their 
corresponding quantal analogues via the correspondence principle. To satisfy 
the reciprocity relation between the cross-section for the direct process 
and o fi for its inverse, we must have 




(19) 


where u)^, and u>^ and u^. are respectively the statistical weights and 
velocities (in a.u.) for the initial and final states of the system. For the 
present problem we have 


w c + 1) *nd s. 


( 20 ) 


From ( 19 ) and (20) we conclude that o^ must have the form 



(Ciwse-sectlonal expression 
symmetric in i and f) 


(21) 
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Extending Kramers type prescriptions for the principal quantum number 
n -* n + 1 to the positive energy continuum where n is replaced by iT], 
Eiedenharn and Brussard (I 965 ) obtain the correspondence 


*\ (n+0 1 


for the so-called Sommerfeld number Tj where 


h X = K \ 2^^ ™ I {23) 

1 , ( fi2 ._£y- 

ViC 

is given in terms of half the distance of closest approach r Q , or the 
incident velocity H . In (25) and z x are the mass and charge of the projectile 
and z 2 stands for the effective charge of the target. From the fundamental 
correspondence between the quantum mechanical matrix elements and the classical 
Fourier components .and from the relation (22), it can be shown (for a concise 
derivation and an elegent discussion see, Biedenharn and Brussard, 1965 ) 
that the classical T| yields the final correspondence 


l J 


note that in terms of the initial and final parameter T)^ and T|^. the quantity 


§ defined in (10) reduces to 


"5 = li-T* 
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For the cross-sectional expression to he inserted in equation (21) we shall 
adopt the classical expression (13 ) and symmetrise it through the use of 
equations (22), (23) and (24) as follows: 

t:, . Un^'XV -')] 4 


.-2A-n. 
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Finally, substituting (26) and (27) in (13) and combining the result with 


(21 ) we obtain 


*5 - S 

|4 1 1 

* 4B A (M c ^n + i f ) Jx(i-«|rlO . 
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An expression similar to (28) is obtained for the differential cross-section 

with f , (§ ) replaced by jjjj, lUi 
A djfl_ 

We shall emphasize the fact that the symmetrization procedure adopted here 
for the classical cross-sections is by no means unique in character. 

Nevertheless it has been found that an essentially similar procedure adopted 
by Biedenham and Brussaard (1965) In nuclear excitation problems, produced 
excellent agreement with the corresponding quantal calculations . We note, 
however, that our symmetrization procedure has the advantage over that of the 
previous authors in that it yields the correct threshold law: ±. Q. 

(Wigner, 1968) for the repulsive Coulomb scattering while the other expression yields 
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. - znh , 

^ ~ — - s £ r near the threshold. The two procedures, however, differ 
little numerically soinewhat away from the threshold. 


4. Numerical Calculations 

The vibrational matrix elements are calculated by using Morse functions 
(Morse, 1929 ) for the target H 2 + molecular ions: 


t«- e--^ k - 2n -'> F„(v 


wher* 


1 Y\ - - 

1 1 \ /L 


vhere x = 


JSa. 

8nDr 


2 = k £ 1 R.® R At 

^ * 

Vtf m iU-,)Q«-0 A 2 m, , ** 

is the spectroscopic parameter, Mq is the v*<iuC*4 mass. 


D q is the dissociation energy and E q1 is the energy difference of the first- 

excited state from the ground state. 

For H 2 + molecular ions we choose 

D q - 0.2053 R 3 M. 

cd # - 0.024ft Ryd. 

M = 918.33^ a.u. 
o 

The energy difference E is calculated by Cohen et al. (i960) and is very 


nearly equal to O.Oe Ryd. 
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Fol lowing Heaps and Herzberg (1952 ) the vibrational matrix elements are 


found to be 


[ t # Oo f 

rnz ) y ^v" t ' y V i ) 2C^ ( V4 TW 
^ n+l 1'- 


where 


q f ^2.V\;-H ) *« It > - fov»i )*»^ ; V 

{ I - t-(n £ .-•■ It-MJiM/ 

and is the polynomial the first few of which are 

H^o = - 0 ! Xn ( I - ("+■* DM a «* 


«•! =0 , a t =i, a 3 = 3 , a 4 -u j o s > JO, 
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5. Results and Discussions 

We present here both the differential as well as the total cross-sections 
for vibro- rotational transitions between the ground state and the first few 
excited vibrational-rotational states. In figures 2 and 3 we compare the 
difference between different rotational excitations and a fixed vibrational 
transition n^= 0 -* n f = 1. We find that the cross-sections for 

= 0 -* = 2 is by far the largest compared to the rest and those between 

= 2, 4, 6 "* jf = 4, 6, 8 respectively are comparable with each other. 

In figures k- and 5 we have similar results between the same set of rotational 
states but with a different vibrational transition: n^ = 0 -* n^ = 2. 

Comparing Fig. 4 and 5 with Fig. 2 and 3 we find that the magnitude of all 
the 0 -* 2 cross s( otions are an order of magnitude or more larger than those 
for the 0 -* 1 transitions. Finally in Fig. 6 we compare the total cross- 
sections for a fixed rotational excitation J^ = 0 -* J^ = 2, with various 
vibrational excitations: n^ = 0 -* n^ = 1, n^ = 0 -* n^. = 2, and 
tt i = 0 -♦ n f = 3 . 

This result again shows that the n^ = 0 -* n^ = 2 vibrational transition 
is much larger than those for n^ = 0 n^ = 1 and n^ * 0 -* n^ = 3j while 
the last two cross-sections are comparable in magnitude. Although we have not 
plotted the cross sections for nj_ = 0 -* n f = 4, we noted that it was somewhat 
smaller than that for n^ = 0 -* nj = 3 in the energy range under consideration. 

We believe that the present method is particularly suitable for the 
energy range E^ = 0.0 to £ 0.5 Ityd. for the following reasons. 

(i) In this energy range the distance of closest approach is ^ 2 a 0 and is 
expected to be large enough to allow essentially non- penetrating orbits 
for the projectiles. 
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(ii) The threshold for positronium formation being at 0.5 Ryd., we may neglect 
such channels fdr positron energies below 0.5 Ryd. 

We conclude by noting that the calculated absolute transition probabilities 
are all very small so that the application of a first order theory is generally 
satisfactory. 
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